On Heatshield Shapes for Mars Entry Capsules 


Dinesh K. Prabhu 1 and David A. Saunders 2 
ERC, Inc., NASA Ames Research Center, Moffett Field, CA, 94035 


The 70° sphere-cone - the standard geometry for all US Mars entry missions - is 
thoroughly examined via flow field simulations at a select few peak heating points along 
candidate flight trajectories. Emphasis is placed on turbulent heating based on the Baldwin- 
Lomax turbulence model. It is shown that increased leeward turbulent heating for a 70° 
sphere-cone flying at angle of attack is primarily due to the discontinuity in curvature 
between the spherical nose cap and the conical frustum - the attachment of the sonic line at 
this sphere-cone junction leads to a supersonic edge Mach number over the leeward acreage. 
In an attempt to mitigate this problem of elevated turbulent heating, alternate geometries, 
without any curvature discontinuities in the acreage, are developed. Two approaches, one 
based on nonlinear optimization with constraints, and one based on the use of non-uniform 
rational B-splines, are considered. All configurations examined remain axisymmetric. The 
aerothermal performance of alternate geometries is shown to be superior to that of the 70° 
sphere-cone. 


Nomenclature 


a = 

Cp 

C Pf 

D b 

L 

M 

Pr 

P 

4 

R 

Re 

R b , R n , R s 

r = 

5c 

St 

T 

V 

x = 

7 

z 

Gi'cek letters 
a = 

r 

A 

s 

£ = 

6, 0 C 

K = 


speed of sound (m/s) 
pressure coefficient 

frozen specific heat at constant pressure (J/kg.K) 
base diameter of heatshield (m) 
characteristic length (m) 

Mach number (= Via) 

frozen Prandtl number (= l uCp f /K) 

pressure (kPa) 

heat flux (W/cm 2 ) 

specific gas constant (J/kg.K) 

Reynolds number (=pVLIp) 

base, nose, and shoulder radii of heatshield (m) 

radial coordinate 

Schmidt number (=/u/pD) 

Stanton number 
temperature (K) 
velocity (m) 
roll axis 

pitch axis, wall-normal coordinate 
yaw axis, radial coordinate 

angle of attack (deg) 

isentropic exponent (or ratio of specific heats) 
shock stand-off distance (m) 
boundary-layer thickness (m) 
emissivity 

angle, cone half angle (deg) 

curvature (m _1 ), thermal conductivity (W/m.K) 


1 Senior Research Scientist, Aero thermo dynamics Branch, 230-3, and Associate Fellow AIAA. 

2 Senior Research Scientist, Aero thermo dynamics Branch, 230-2, and Senior Member AIAA. 

1 

American Institute of Aeronautics and Astronautics 



p 

p 

r = 

co = 

Subscripts 
00 = 

c = 

e = 

t 

w = 

Superscripts 
* = 


viscosity (Pa.s) 
mass density (kg/m 3 ) 
shear stress (Pa) 
vorticity (s' 1 ) 

freestream 

compressible 

boundary-layer edge 

turbulent 

wall 

non-dimensional quantity 


I. Introduction 

T HE success of NASA’s Viking mission in the 1970s has meant that all subsequent missions to Mars have relied 
on a 70° sphere-cone capsule configuration to land scientific experiments on the surface of the planet; staying 
with a proven architecture/configuration carries less risk. The paper of Braun and Manning 1 provides a good 
summary of all NASA missions to Mars starting with Viking, and the challenges of entry, descent and landing 
(EDL) at Mars. Table 1 lists the entry type and heatshield geometries (along with key parameters) of capsules used 
by the United States (US) in several Mars missions from 1975 through 2011. The data have been extracted from the 
compilation of Davies and Arcadi. 2 


Table 1. Heatshield parameters for several US Mars missions. 



Entry Type 

HS Geometry 

D b 

m 

Rn 

m 

R s 

m 

R,/R„ 

R/R„ 

Deep Space 2 a 

Ballistic 

45° Sphere-cone 

0.35 

0.088 


0.500 


Viking 

Lifting 

70° Sphere-cone 

3.54 

0.893 

0.025 

0.500 

0.014 

Mars Pathfinder 

Ballistic 

70° Sphere-cone 

2.65 

0.663 


0.500 


Mars Exploration Rovers 

Ballistic 

70° Sphere-cone 

2.65 

0.663 


0.500 


Mars Phoenix 

Ballistic 

70° Sphere-cone 

2.65 

0.663 

0.066 

0.500 

0.050 

Mars Science Laboratory 6 

Lifting 

70° Sphere-cone 

4.50 

1.125 

0.125 

0.500 

0.056 


a Mission presumed failed 

b During writing of the present paper, the MSL vehicle was launched successfully on Nov. 26, 2011 


The latest in the series of Mars missions is the Mars Science Laboratory (MSL), which was launched on 
November 26, 2011. 3 This 4.5 m diameter 70° sphere-cone configuration will be NASA’s largest capsule to enter 
the Mars atmosphere (or any planetary atmosphere for that matter). Like Viking which flew a lifting trajectory with 
a nominal angle of attack of -11°, the MSL capsule too will fly a lifting trajectory, only this time with a nominal 
angle of attack of -16°. 

As with any blunt capsule at angle of attack, larger size means longer running lengths of streamlines (from the 
stagnation point) towards the leeward shoulder and hence increased likelihood of transition to turbulence, with a 
consequent increase in heating on the leeward side. 4 For certain regions of the flight trajectory space (parameterized 
by entry velocity and entry flight path angle), MSL is likely to experience significant turbulent heating, a situation 
that will be further exacerbated by surface roughness that develops due to ablation of the Thermal Protection System 
(TPS). To illustrate the nature of the problem, contours of computed (and unmargined) surface turbulent heat flux, 
along with centerline distributions at the peak heating point along a flight trajectory for a 4.7m diameter 70° sphere- 
cone (a candidate geometry for a future Mars mission) are shown in Fig. la and Fig. lb. On the leeward side, the 
elevated levels of turbulent heating, as much as 6 to 7 times the corresponding laminar value, are clearly seen. In the 
case of MSL, this likelihood of transition and heating augmentation due to surface roughness, coupled with material 
failure observed during arc -jet testing, 5 led the project to abandon SLA-561V (Super Light Ablator) as the first 
choice of thermal protection material with Mars Pathfinder flight heritage, and use PICA (Phenolic-Impregnated 
Ceramic Ablator) instead. 6 

Despite mitigating the risk associated with turbulent heating, one question that has been left unanswered in the 
MSL project is: are there other heatshield shapes that result in lower turbulent heating on the leeside? There is 
currently a debate within the US Mars program community on a change to the heatshield shape. The first choice for 
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a new heatshield for Mars entry is naturally one with flight heritage (albeit for Earth entry), and that is a spherical 
section heatshield used in several space programs including Apollo and Orion. 

In an attempt to answer the question yet unanswered by the Mars program, a fresh look at heatshield shapes has 
been taken in the present work, and a few “non-traditional” (yet analytic) shapes have been developed and analyzed 
using modem flow field simulation tools. 7 8 It should be noted that the idea of moving away from “traditional” 
capsule shapes - sphere-cones or spherical sections - is not new. For instance, as part of NASA’s Exploration 
Systems Architecture Study (ESAS), 9 asymmetric capsule shapes were developed and studied. 10 Although the 
“traditional” symmetric spherical section heatshield (Apollo shape) was finally selected for the Orion capsule, the 
groundwork had been laid for systematic development and analysis of asymmetric configurations defined by a finite 
number of geometric parameters. The present paper builds upon this foundation, and new symmetric shapes that lie 
between a sphere-cone and a completely spherical section are developed and studied. Retaining axisymmetry still 
leaves plenty of room for variation, so any possible benefits of asymmetry are not pursued here. The focus of the 
present computational study is only on the aerothermodynamic environments for these new shapes. It is hoped that 
these configurations, after more experimental testing, will be considered as possible choices for future MSL-class 
Mars missions, or perhaps as starting points for more elaborate shape optimization studies. 
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(a) Surface heat flux (b) Centerline heat flux 

Figure 1. Surface and centerline distributions of hot-wall ( e = 0.85) heat flux for a 70° sphere-cone at -15.7° 
angle of attack for a freestream velocity of 5.41 km/s. Elevated levels of turbulent heating (Baldwin-Lomax 
model) on the leeward side are clearly seen. Results are for a smooth wall, and the wall is assumed 
“supercatalytic” to recombination. 


II. Ground Rules and Assumptions 


A. Geometries 

As already indicated, only axisymmetric forebody configurations are considered in the present work. This choice 
simplifies the construction of three-dimensional candidate heatshield geometries. One only needs to construct the 
generatrix (a 2D plane curve) and simply rotate that generatrix around the x-axis (considered the capsule roll axis in 
the present work) to construct the three-dimensional surface representing the outer mold line (OML) of the 
heatshield. The two most commonly used geometries for axisymmetric heatshields of atmospheric entry capsules are 
the sphere-cone and the spherical section whose generatrices are respectively shown in Figs. 2a and 2b. 

The generatrix of a sphere-cone heatshield (Fig. 2a) is actually composed of three segments - a circle segment 
for the nose, a straight line segment for the cone, and another circle segment for the shoulder - which are blended 
smoothly, i.e., there is no discontinuity in slope where any two segments meet. This composite plane curve is 
completely described by four parameters - R b , R m R s , and 6 C - where R n < R b , R s « R b (or R n ), and 0° <6 C <90° 
(70° is usually the limit!). 
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The generatrix of a spherical section heatshield is composed of two segments (Fig. 2b) - a circle segment 
representing the dish, and a circle segment representing the shoulder. Again, the two segments are blended smoothly 
at the dish-shoulder junction. This second composite plane curve is completely described by three parameters - R b , 
R m and R s - where R n > R b , and R s « R b (or R n ). 



(a) Sphere-cone (b) Spherical section 

Figure 2. Generatrices of typical symmetric heatshield geometries: (a) a sphere-cone parameterized by the 
radius, R n (R n < R h ), of the nosecap, the radius, R s (R s « R b or R lt ), of the shoulder torus, and the cone half 
angle, 6 C ("0° < 0 C < 90°); and (b) a spherical section parameterized by the radius, R n (R n > R b ), of the dish, and 
the radius, R s (R s « R b or R n ), of the shoulder torus. The base radius for either geometry is R b . 


Both the composite curves shown in Figs. 2a and 2b are simple enough that analytical expressions for the 
coordinates of points of intersection of the component segments can be obtained with straightforward algebra. For 
instance, if the apex of the sphere-cone heatshield generatrix (Fig. 2a) is at the origin (0,0) of the coordinate system 
(x,r), then the coordinates (x t>m r t>n ) of the point of tangency between the circle and the straight line (inclined at 6 C to 
the x-axis) are: 


X Un =^,[ 1 - Sin6 c] r ,,n = R n COS0 c ( la ) 

and the coordinates (x t>s , r t>s ) of the point of tangency of the straight line and the shoulder circle are: 

=^(l-sin0 c ) + [^ -R n COS0 C -fl s (l-cos0 c )]cot0 c r ts = R b -/?,[l-cos0 c ] (lb) 

Similarly, for the spherical section generatrix (Fig. 2b), the coordinates (x tiS , r t>s ) of the tangency point between 
the circle segment representing the dish and the circle segment representing the shoulder are: 


V, = R n 


1 - 


(R n -R b ){R n+ R b -2R) 


R, - R „ 


V 


r = * ^ 

us " R„ - R. 


(2a) 


The angle, 6, between the x-axis and the tangent to the dish at the intersection point with the shoulder is: 


6 = arccos 


R „- R s 

R„ -R. 


(2b) 


An important geometric aspect of a heatshield, and one often overlooked in flowfield analysis, is the distribution 
of curvature of the three-dimensional surface. The choice of axisymmetric shapes considerably simplifies the 
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determination of curvature. One needs to consider only the radial distribution of the one-dimensional curvature of 
the generatrix (a plane curve), i.e., the longitudinal curvature of the three-dimensional surface swept by the 
generatrix; the other principal curvature of the surface, the one in the circumferential direction, is a constant for a 
given x value due to the choice of axially symmetric shapes. The radial distribution of one-dimensional curvature 
(magnitude) of the generatrix is: 






3 

K(r) = 

d 2 x 


' dxY 

~2 

dr 2 

1 + 

\ dr ) 








The distributions of curvature for the two most common generatrices, the sphere-cone and the spherical section 
(Figs. 2a and 2b), are shown in Fig. 3. 



Figure 3. Curvature ( k ) distributions of sphere-cone and spherical section geometries. The abbreviation SC 
stands for sphere-cone, and SS stands for spherical section. For either geometry, curvature changes 
discontinuously at the points of tangency between component segments. 

It should be noted that while the geometric construction (Eq. 1 or Eq. 2) guarantees continuity of the first 
derivative in r between various segments of the heatshield, it does not guarantee continuity in the second derivative, 
and hence, curvature. It is hypothesized that a discontinuity in the second derivative would have an impact on 
predictions of surface aero thermal environments for these traditional shapes. This hypothesis is put to the test in the 
present paper. 

From Table 1 it is evident that, with the sole exception of Deep Space 2, all US Mars missions with an 
atmospheric entry component (lifting or ballistic) in the mission profile have used a 70° sphere-cone geometry for 
the heatshield. The 70° sphere-cone shape has the high drag necessary to decelerate the capsule through the tenuous 
Martian atmosphere. Therefore, the point of departure for a heatshield shape study - the focus of the present paper - 
is a 70° sphere-cone. 

For all the cases listed in Table 1, the radius of the spherical nose cap is always one half the base radius, i.e ., 
RJRb = 0.5, and the shoulder radius is roughly one tenth that of the nose, i.e., RJR n ~ 0.1 or RJRb ~ 0.05. There is no 
information in the extant literature for the rationale behind the particular values of these ratios, especially RJRb . The 
influence of these ratios is also explored in the present work. 

Although a spherical section heatshield has never been used in US Mars missions, such a shape is the baseline 
for human-rated Earth-entry capsules such as Mercury, Gemini, Apollo, and the current Multi-Purpose Crew 
Vehicle (MPCV). Therefore, the spherical section heatshield is considered as an alternate baseline in the present 
work. The aerothermal performance of other symmetric shapes developed in the following sections is gauged 
against these two shapes with flight heritage. 
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B. Modeling and simulation 

The modeling assumptions used in the present study are identical to those used in the development of 
aerothermal design environments for MSL. 3 Key modeling assumptions used in LAURA 6 and Dplr 1 simulations for 
MSL are summarized in Table 2. In the present work version 4.02.2 Dplr is used in all simulations; v3.05 of Dplr 
was used in MSL simulations. 


Table 2. Models used in flow field simulations with v4.02.2 of Dplr. 


Model 


Details 


Gas mixture 


Thermodynamics 


Laminar transport 


Turbulent transport 


Reaction mechanism & rates 


Wall boundary condition 


8-species (C0 2 , CO, N 2 , 0 2 , NO, C, N, and O) Mars gas mixture 
A 97-3% mixture, by mass, of C0 2 -N 2 is assumed as the freestream 
Two-temperature ( T-T v ) model 

Simple Harmonic Oscillator/Rigid Rotor (SHO/RR) for specific heat and enthalpy of 
each molecular species 

Excited electronic states are also included for all species (atoms and molecules) 

Collision-integrals for species pairs with Gupta-Yos mixing rule 

Mass diffusion modeled by SCEBD (Self-Consistent Effective Binary Diffusion) 

• The SCEBD model, an approximation to a true multi-component diffusion model, 
does not require specification of a Schmidt number 

Baldwin-Lomax algebraic model including modifications for compressibility and 

stream wise pressure gradient 

Prandtl number ( Pr t ) = 0.9 

Schmidt number (Sc t ) = 0.7 

Mitcheltree mechanism 

• Rates from Park et al. 

“Supercatalytic” to recombination 

• “Supercatalytic” means that the freestream composition is recovered at the wall, i. e . , 
Dirichlet boundary conditions are used for mass fractions at the wall 

Radiation equilibrium with e = 0.85 

• Energy balance between heat brought to the wall by convection and diffusion and re- 
radiated into free space (at 0 K). No heat conduction through the thickness of the 
material 


The turbulence modeling aspect requires a little further elaboration. The turbulence model used in Dplr 
simulations for MSL and other geometries discussed here is the zero-equation algebraic model of Baldwin and 
Lomax. 11 This model, originally developed for incompressible flows, has been modified to account for flow 
compressibility 12 and the influence of adverse pressure gradients in the stream direction. 13 These two modifications 
to the inner layer of the Baldwin-Lomax model are: 


Vc =r l 




A + =26 


(4) 


where rj is the directed distance from the wall, r w is the shear stress at the wall, p t is the turbulent eddy viscosity, 
and \co\ is the magnitude of vorticity. It should be noted that the turbulent eddy viscosity enters into the model 
implicitly in the modified definition of A + . 

The influence of compressibility is accounted for when local values (instead of wall values) of p and p are used 
in the definition of rf . The influence of the streamwise pressure gradient comes in through the parameter A + , which 
usually has a value of 26. It is scaled by the ratio of the shear stress, r W9 at the wall to the local value, r. The 
modified rf and A + are then used in the Van Driest damping function, D\ 


( „ + YI 


D = 


1 - exp 


^7c 

a: 


l mix= K VD K = 0.4 


P J . 


(5) 


Although the correction for compressibility has an influence on the magnitude of D in Eq. 5, the influence of the 
pressure gradient correction term is larger, especially since it includes the magnitude of local vorticity and the local 
turbulent eddy viscosity. Any decrease in the exponential term means increased mixing length, and hence increased 
turbulent eddy viscosity because the turbulent eddy viscosity depends on the square of the mixing length. 
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Furthermore, depending on the curvature of the shock, its standoff from the body, and the Reynolds number (which 
dictates the thickness of the wall-bounded shear layer), the attached shear layer can entrain rotational flow generated 
by a curved shock (Crocco’s theorem 14 ), and hence cause enhanced heating. A second line of inquiry in the present 
paper is the interplay between the turbulence model and discontinuities in surface curvature. 

C. Flow conditions for flight simulations 

It is well known that the maximum predicted heat flux, with appropriate margins to account for uncertainties in 
modeling and simulation, along a flight trajectory dictates the choice of TPS material for the heatshield. 15 Therefore, 
only peak heating points along a few design trajectories for a proposed 2018 US mission to Mars are considered in 
the present work. The freestream conditions necessary for CFD simulations at these points are given in Table 3. The 
thickness of the chosen TPS material, which is dependent on the time-integrated value of heat flux (i.e., heat load) 
over a trajectory, is not considered here. 


Table 3. Freestream conditions chosen for present study. 


Trajectory 

V x 

m/s 

a 

deg 

P OO 

kg/m 3 

xlO 4 

T x 

K 

M x 

Re «, 
1/m 
xlO 6 

Qc O 

kPa 

a 

q oo 

W/cm 2 

Retrograde/Steep 

5411.2 

-15.71 

10.13 

182.6 

24.2 

0.585 

14.83 

96 

Prograde/Nominal 

5070.9 

-15.54 

9.06 

175.9 

23.3 

0.510 

11.65 

75 

Prograde/Shallow 

5165.0 

-15.54 

6.87 

168.8 

24.2 

0.411 

9.16 

69 


a Cold-wall heat flux, computed using Sutton-Graves formula [Ref. 16], at the stagnation point of a hemisphere of 1 m radius 


III. Preliminary computations with traditional shapes 

First, baseline heating environments are established for the two candidate geometries - a 70° sphere-cone 
geometry and a spherical section geometry. Both heatshield geometries have a base diameter of 4.7 m, a nearly 5% 
increase over the diameter of MSL. At the time when the present study was undertaken, mission planners had 
considered an increase over the diameter of the MSL geometry to accommodate a larger payload. 

Next, since the focus of the present paper is on turbulent heating over the heatshield, we calibrate the simulation 
tool - v4.02.2 of Dplr - to ground-based turbulent flow experiments (with and without any real-gas effects). Having 
calibrated the simulation tool, we then systematically examine the influence of real-gas effects, bluntness, curvature 
continuity (or lack thereof), angle of attack, freestream conditions, etc., on turbulent heating before developing 
alternate geometries for improved aerothermo dynamic performance of the heatshield under turbulent flow 
conditions. 

A. Establishing the baseline: 70° sphere-cone vs spherical section 

For the peak heating point along a retrograde/steep entry trajectory (Table 5), laminar and turbulent heating 
computations are performed for both the 70° sphere-cone and spherical section geometries. The centerline 
distributions of hot-wall (£=0.85) heat flux and pressure for the two geometries are shown in Fig. 4. 
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(a) Heat flux (b) Pressure 

Figure 4. Centerline distributions of hot-wall (s = 0.85) heat flux and pressure for a 70° sphere-cone and a 
spherical section heatshield (both 4.7 m in diameter). Laminar and turbulent (Baldwin-Lomax) results are 
shown for a flight speed of 5.41 km/s and an angle of attack of -15.7°. Dashed lines represent laminar results. 

For the laminar case, the spherical section has heating levels that are roughly no more than 15% of the levels on 
the windward and leeward flanks of the cone (Fig. 4a). However, the spherical nose of the cone experiences much 
larger heating. For the turbulent case, the spherical section experiences much less heating (about 50% less) over 
much of the leeward acreage. Unlike the spherical section, which has a smooth distribution of turbulent heating, the 
sphere-cone experiences a very rapid rise in turbulent heating level past the sphere-cone junction on the leeward 
side. 

A probable cause for this lies in the pressure distributions shown in Fig. 4b. The spherical section has a favorable 
pressure gradient over the entire wetted acreage, while the sphere-cone pressure distribution is influenced by the 
curvature discontinuity at the junction of the spherical nose cap and the conical frustum. Furthermore, the laminar 
and turbulent pressure distributions towards the leeward side differ significantly in character for the sphere-cone 
geometry compared to the spherical section, which does not show a difference in pressure distributions for laminar 
and turbulent flow. For the 70° sphere-cone geometry, the laminar pressure distribution beyond the point of 
curvature inflection is smoother on the leeward side, but the pressure gradient appears to be somewhat unfavorable 
in the streamwise direction. The turbulent pressure distribution, on the other hand, has multiple inflection points 
towards the leeward shoulder. These inflection points appear to line up with rate of change of heat flux towards the 
leeward shoulder. This behavior of surface pressure is of great interest, and is investigated systematically in the 
present section. 

It is first necessary to examine available experimental data obtained from various ground-based facilities to see if 
the measured and predicted heating distributions correlate with the predicted pressure distributions and inflections 
(if any) in them. For the 70° sphere-cone geometry, there is a wealth of aerothermal testing data acquired under the 
MSL program. There is also a wealth of test data acquired for the spherical section under the MPCV (Multi-Purpose 
Crew Vehicle) program. However, the spherical section geometry developed under the Mars 2018 program differs 
from that of the MPCV in terms of the dish radius, R m and the shoulder radius, R s . Therefore, there are no 
experimental data for that particular shape. 

B. Ground-based experiments 

The first test case, that of high Reynolds number flow of pure N 2 , is one taken from a matrix of conditions 
presented in the comprehensive review paper of Hollis. The freestream conditions for this case are presented in 
Table 4. At these conditions, a 15.24 cm diameter 70° sphere-cone (0.0339-scale MSL) and a 17.78 cm diameter 
spherical section (0.0354-scale MPCV) were tested in the T9 tunnel (a conventional blowdown tunnel) at Arnold 
Engineering Development Center (AEDC) in Silver Spring, MD. The angles of attack for the two shapes, however, 
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were not identical in the tests - the 70° sphere-cone was tested at an angle of attack of -16°, while the spherical 
section was tested at an angle of attack of -28°. Although the freestream Reynolds number is high enough to 
guarantee fully turbulent flow over the entire heatshield (see Ref. 17) for both shapes, the stagnation enthalpy of the 
flow is quite low. Therefore, the simulations of this case address only the fluid mechanical aspects of the turbulent 
flow problem and are not encumbered by any real-gas effects. 


Table 4. Freestream condition for a case from the AEDC T9 Tunnel MSL/MPCV test matrix. 


Test 

M „ 

/?e m xl0 7 

Poo 

Tc „ 

Rco 

t/oo 

H o-Ah^ook 

Gas 


(1/m) 

(kPa) 

(K) 

(kg/m 3 ) 

(m/s) 

(MJ/kg) 

n 2 

8.0 

16 

11.8553 

70.1 

0.5718 

1356.8 

0.681 


For the freestream condition given in Table 4, flow computations are performed, first assuming the flow to be 
completely laminar, and then assuming the flow to be fully turbulent. In all computations, the wall is assumed 
isothermal with a temperature of 300 K. Since the flow is of pure N 2 and the stagnation enthalpy, H 0 , is low, 
computations are performed with y = 1.4 and a specific gas constant, R , of 296.78 J/kg.K. Sutherland’s law is used 
for molecular viscosity, and the thermal conductivity of the gas is computed with the assumption of a constant 
Prandtl number (Pr = 0.71). The Baldwin-Lomax algebraic turbulence model, corrected for compressibility effects 
and streamwise pressure gradient, is used in the turbulent flow simulations. It is assumed that the flow transitions 
from laminar to turbulent at the stagnation point, and therefore, no model for transition is employed in the turbulent 
flow simulations. 

The predicted centerline distributions of cold-wall heat flux (laminar and turbulent) are compared against 
experimental data in Figs. 5a and 5b for the two shapes, 70° sphere-cone and spherical section, respectively. Both 
predictions and experimental data are nondimensionalized as: 



Figure 5. Centerline non-dimensional cold-wall (T W =300K) heat flux (or equivalently, scaled Stanton number) 
distributions for: (a) a 70° sphere-cone at -16° angle of attack, and (b) a spherical section at -28° angle of 
attack. Both configurations were tested at a Mach number of 8 in high Reynolds number flow of pure N 2 in 
the T9 tunnel of AEDC. 


The turbulent flow predictions using the Baldwin-Lomax algebraic turbulence model are in good agreement with 
the experimental measurements. As reported in the work of Hollis, 17 the experimental data have an uncertainty of 
±12%. For the most part, predictions lie within this uncertainty band. It must be noted that computations for both 
shapes have not been performed with excessively fine grids. However, meshes in either case have been tailored to 
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align with the bow shocks and employ sufficient grid point clustering in the wall-normal direction to ensure good 
resolution of the wall-bounded shear layers. 

Predicted centerline distributions of pressure are shown in Figs. 6a and 6b, for 70° sphere-cone and spherical 
section, respectively. 




Figure 6. Centerline pressure distributions for: (a) a 70° sphere-cone at -16° angle of attack, and (b) a 
spherical section at -28° angle of attack. Both configurations were tested at a Mach number of 8 in high 
Reynolds number flow of pure N 2 in the T9 tunnel of AEDC. 


In the case of the 70° sphere-cone, the discontinuity in curvature (see Section IIA) at the sphere-cone junction 
leads to an inflection in the pressure distribution, regardless of whether the flow is laminar or turbulent. Beyond this 
junction, the streamwise pressure gradient is favorable ( i.e ., the gradient has a negative sign) towards both the 
windward and leeward shoulders. In contrast to this behavior, the spherical section, which does not have any 
curvature inflection apart from at the shoulder, has a favorable pressure gradient over the entire leeward side of the 
heatshield. This distribution is similar to the distribution at -15.7° for the flight case (Fig. 4b), but steeper because of 
the higher angle of attack (-28° vs -15.7°). For either configuration, the pressure distributions predicted for laminar 
and turbulent flows are identical. 

Since the T9 test case does not have any high-enthalpy effects (or real-gas effects), a second test case, that of 
high Reynolds number flow of dissociated C0 2 , is one taken from a matrix of conditions presented in the paper of 
Wright et al. n The freestream condition for this case are presented in Table 5. At this condition, a 17.78 cm 
diameter 70° sphere-cone (0.0395-scale MSL) model made of stainless steel was tested in the T5 tunnel (a reflected 
shock tunnel) at Caltech. The angle of attack of the model was -16°. For the freestream condition given in Table 5, 
the freestream Mach and unit Reynolds numbers are M ^ = 4.61, Re ^ = 0.877 xlO 7 1/m, respectively. The unit 
Reynolds number for this case is about 20 times lower than in the T9 case. 


Table 5. Freestream condition for a case (Shot #2263) from the Caltech T5 Tunnel MSL test matrix. 


Ho 

(MJ/kg) 

Po 

(MPa) 

Coo 

(m/s) 

Poo 

(kg/m 3 ) 

Too 

(K) 

c[C0 2 ] 

c[0 2 ] 

c[CO] 

C[ O] 

5.7 

71 

2701 

0.17769 

1470 

0.898 

0.037 

0.065 

0.0 


This second case differs from the first in two ways. Firstly, the freestream is a high-temperature dissociated gas 
mixture consisting of C0 2 , CO, and 0 2 . This mixture is further heated by the bow shock that envelopes the test 
article. As a consequence, y is no longer a constant, and is a function of both temperature and chemical composition. 
Secondly, chemical reactions at the surface of the test article are important, even if the wall is assumed isothermal at 
300 K. Chemical reactions at the surface are recombination of atomic oxygen diffusing down to a cold surface, and 
oxidation of CO. The latter is a surface-mediated reaction and often requires a noble metal such as platinum or 
palladium to catalyze it. 
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A simple approach towards surface chemistry has been adopted in the development of environments for 
aerothermodynamic design (TPS material selection and sizing) of the heatshield of the MSL. The wall is assumed to 
be “supercatalytic” to recombination of atomic and molecular species, 19 i.e., the freestream composition is 
completely recovered at the surface, and hence all the enthalpy in the chemical modes of the dissociated gas 
mixture. Since this “supercatalytic” model provides conservative estimates of surface heating, it is used for design 
purposes. For the case of a dissociated freestream, the term “supercatalycity” becomes ambiguous because recovery 
of the freestream composition at the body surface would mean that chemical enthalpy recovery is incomplete. In the 
work of Wright et al . 18 this ambiguity was removed by assuming that the chemical composition at the body surface 
was 100% C0 2 regardless of the state of dissociation of the freestream. In the present work, the actual freestream 
composition is recovered at the surface, and it should be expected that predicted heat fluxes for this ground test 
would be lower than those in the work of Wright et al. 

An alternate approach to the problem of surface chemistry is to assume that there are two competing rate 
processes for atomic oxygen adsorbed on the surface 20 - one for recombination with another oxygen atom in the gas 
phase, and one for recombination with CO in the gas phase, i.e., 

(1) O + 0[s] — > 0 2 

(2) C0 + 0[s]^C0 2 

This model, dubbed “fully catalytic,” has had some success in replicating heating measurements made in the 
ballistic range. 21 Heat fluxes inferred from thermal imaging of a projectile (a sphere) flying through quiescent C0 2 
(in the ballistic range) are in fair agreement with predictions based on the rate chemistry model. In the present work, 
both approaches have been used in numerical simulations of the test case, and the results are shown in Fig. 7. 




Figure 7. Centerline distributions of: (a) non-dimensional cold-wall (7k,=300K) heat flux, and (b) non- 
dimensional pressure for a 17.78 cm diameter 70° sphere-cone tested at a Mach number of 4.6 and at an angle 
of attack of -16° in high Reynolds flow of dissociated C0 2 in the reflected shock tunnel T5 at Caltech. 

Predicted heat fluxes along the centerline of the test model are compared against measurements in Fig. 7a. An 
arbitrary ±15% error bar is used to indicate uncertainty in measurement. The heat fluxes predicted by the 
“supercatalytic” wall model are in good agreement with measured data, while those predicted by the “fully catalytic” 
wall model are significantly lower. However, the trends in turbulent heating on the leeward flank of the cone are 
similar for both predictions and measurements. Although not shown in Fig. 7a, the heat fluxes predicted by the 
“supercatalytic” wall model employed in the present work are lower than those predicted by Wright et al. (see Fig. 
15 in Ref. 18). This is as expected because of the ambiguity in the meaning of the term “supercatalytic”. The 
variation of pressure along the centerline of the model is shown in Fig. 7b. 

Results of numerical simulation of ground-based experiments in two different facilities make for interesting 
comparisons. The distributions of heat flux for the 70° sphere-cone tested in T9 tunnel at AEDC (Fig. 5a) and the 
reflected shock tunnel in Caltech (Fig. 7a) are quite similar in that the leeward turbulent heating levels are nearly 
“uniform.” However, the pressure distributions (Figs. 6a and 7b) show differences. Firstly, the pressure expansion at 
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the sphere-cone junction is shallower for the T9 case compared to the T5 case. Secondly, while the stream wise 
pressure gradient is favorable toward the leeside for the T9 case, it is not favorable and neutral at best for the T5 
case. Inflections in the pressure distribution and the sign of the pressure gradient are important aspects, and 
discussions about them are deferred until other aspects of the problems, viz., real-gas effects, angle of attack, 
bluntness, etc. are studied in detail. 

Contours of Mach number in the pitch plane of the 70° sphere-cone models tested in the T9 and T5 tunnels are 
shown in Fig. 8. Only three contour levels - 0.01, 1.0, and maximum freestream Mach number - are shown for each 
case. 



(a) T9 Model (15.24 cm diameter), = 7.95 (b) T5 Model (17.78 cm diameter), M x = 4.6 

Figure 8. Contours of Mach number in the pitch plane of the 70° sphere-cone models tested in the T9 and T5 
tunnels. Only three contour levels - 0.01, 1.0, and maximum freestream Mach number (7.95 or 4.6) - are 
shown. Both models are at -16° angle of attack. 

For both the T9 and T5 models, the sonic line attaches at the shoulders (windward and leeward). However, for 
the T5 case, the sonic line has an inflection in the vicinity of the apex because the shock is drawn in closer to the 
body because of the compression of the test gas (C0 2 ). The behavior of the sonic line in the T5 case is still different 
from the flight case (Fig. lc), where due to further compression, the sonic line attaches on the sphere-cone interface 
on the leeward side. 

The pressure distributions from laminar and turbulent flow computations are identical to each other for both the 
T9 and T5 cases. This is unlike the distributions for the flight case (Fig. 4b). Could this behavior be attributed to 
differences in real-gas effects (catchall for excitation of internal structure of atoms and molecules and/or chemical 
reactions) between ground and flight? If real-gas effects are indeed the cause of the differences, then the issue of 
ground-to-flight traceability becomes a little more difficult to address satisfactorily. 

C. 70° sphere-cone: Real-gas effects 

In order to study the influence of internal structure (electronic or rovibronic) of constituent atoms and molecules, 
and chemical reactions between them, a quick computational experiment is performed at the peak heating point of a 
retrograde/steep entry trajectory (see Table 3). Calculations are performed for three variations in the 
thermochemistry model. These models are 
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Model 1 : No dissociation is allowed for the C0 2 -N 2 mixture. Further, internal (vibronic) structure of C0 2 and N 2 
are not considered. Since C0 2 is a linear molecule, disallowing any vibration means that the gas mixture behaves as 
an ideal one with 7 = 1 . 4 , but with a different specific gas constant. 

Model 2 : Dissociation is still not allowed, but the vibrational modes of C0 2 and N 2 are included in the 
thermodynamic model. In this case, y is no longer a constant and depends on the temperature of the gas mixture. 
Thermal equilibrium is assumed, and GRC curve fits 22 are employed for thermodynamic properties. 

Model 3 : The C0 2 -N 2 mixture is allowed to dissociate behind the shock and internal structure of the constituent 
atoms and molecules are considered. For this particular model, the wall is assumed to be noncatalytic to 
recombination in order to keep the boundary condition consistent with the isothermal wall boundary condition used 
for the other 2 models. Additionally, thermal nonequilibrium ( T ± T v ) is assumed in this model. 

Figure 9 shows the centerline distributions of cold- wall (T w = 400 K) heat flux and pressure obtained from 
laminar and turbulent flow computations for each of these three thermochemistry models. 

Lee Wind Lee Wind 




Figure 9. Centerline distributions of: (a) cold-wall (7k,=400K) heat flux, and (b) pressure for a 4.7 m diameter 
70° sphere-cone at a flight Mach number of 24.4 and at an angle of attack of -15.7°. The freestream is a 97- 
3% (by mass) mixture of C0 2 -N 2 . Three variations in the thermochemistry model are considered: (i) without 
internal structure and without reactions, (ii) with internal structure and without reactions, and (iii) with both 
internal structure and chemical reactions. Dashed lines represent laminar results. 

The distributions of laminar heat flux are more or less similar in trend for the three thermochemistry models, 
even if the magnitudes are different (Fig. 9a). However, the distributions of turbulent heat flux differ significantly 
depending on the thermochemistry model employed. An examination of the pressure distributions (Fig. 9b) makes 
clear the influence of internal structure and chemistry. The behavior of pressure beyond the sphere-cone junction is 
different for each of the three thermochemistry models that have been employed. For Model 1 ( 7 = constant), there is 
a small inflection in pressure, but the pressure recovers beyond the expansion and decreases towards either shoulder, 
i.e., the pressure gradient is favorable. Furthermore, the pressure distributions for the laminar and turbulent cases are 
indistinguishable for this model. For Model 2 (7 ^ constant, but no chemical reactions), the depth of the pressure 
inflection at the sphere-cone junction is larger than in Model 1. The pressure distributions for the laminar and 
turbulent cases, while a little different from each other, do not show secondary inflections with large amplitudes. 
Consequently, the rise in turbulent heating is not as dramatic. With Model 3 ( 77 ^ constant and chemical reactions are 
allowed to occur), the pressure distributions for laminar and turbulent flow differ, and inflections in pressure past the 
sphere-cone junction have larger amplitudes. 

We expect that with each refinement in the model starting with Model 1, we allow the shock layer to compress 
more, hence drawing in the bow shock closer to the body. Contours of Mach number in the pitch plane are shown in 
Fig. 10. 
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Axial coordinate, m 


Figure 10. Pitch plane Mach number (24.4 and 1.0) contours for a 4.7 m diameter 70° sphere-cone at a 
flight Mach number of 24.4 and at an angle of attack of -15.7°. The freestream is a 97-3% (by mass) mixture 
of C0 2 -N 2 . Three variations in the thermochemistry model are considered: (i) without internal structure and 
without reactions, (ii) with internal structure and without reactions, and (iii) with both internal structure and 
chemical reactions. 

As expected, the bow shock is drawn in closer to the body by changing the thermochemistry from an ideal-gas 
model (y = constant) to a real-gas model. The sonic lines for each of these models attach at nearly the same point on 
the windward shoulder. However, for the y = constant case, the sonic line attaches on the leeward shoulder, while it 
attaches closer to the sphere-cone junction for the other two cases. Furthermore, the shock standoff distance 
decreases on the leeward side causing the bow shock to acquire a larger curvature in the vicinity of the apex of the 
sphere-cone. The more curved the shock is, the more likely it is to generate vorticity (Crocco’s theorem 14 ) and this 
shock-generated vorticity will likely get entrained within the wall-bounded shear layer - the notion of “entropy 
swallowing.” Vorticity entrainment (or “entropy swallowing”) is an issue for boundary-layer methods because 
specification of edge conditions becomes difficult in the presence of such entrainment. However, this is not an issue 
for Navier-Stokes equations, which are uniformly valid in the shock layer. Nevertheless, we take cognizance of the 
possibility of rotational flow entrainment within the wall-bounded shear layer. 

We hypothesize that the behavior in pressure depends strongly on the shock- standoff distance and attachment of 
the sonic line on the leeward flank of the cone. A reduced shock standoff distance in the vicinity of the apex of the 
sphere-cone would mean increased shock curvature (the flow “feels” the presence of the nose cap), and hence 
generation of vorticity that will likely get entrained within the wall-bounded shear layer. A way to check this 
hypothesis is to blunt the cone, i.e., increase the radius, R n , of the spherical nose cap. By blunting the cone, the shock 
standoff is increased, and the likelihood of vorticity getting entrained within the wall-bounded shear layer is 
reduced. In addition, the attachment of the sonic line on the leeward side of the cone can be studied as well. 

D. 70° sphere-cone: Influence of blunting 

As mentioned earlier, a rule of thumb used in the design of sphere-cone heatshields is to let the nose radius be 
one half the base radius, i.e., RJRb = 0.5. There does not seem to be any evidence in the extant literature that this 
ratio is based on any theory, but more so on the success of missions that have employed capsules based on this rule. 
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Consequently, we relax this assumption and consider blunter noses for which RJRb = 0.8 and 1.1. The generatrices 
for these shapes are shown in Fig. 11. Also shown in the figure is the generatrix for the spherical section heatshield. 



Axial coordinate (x), m 


Figure 11. Generatrices for a 70° sphere-cone heatshield with various values of nose bluntness, i.e., RJRb 
values. Also shown in the figure is the generatrix for a spherical-section heatshield. All profiles are lined up at 
the plane of maximum diameter. 



Figure 12. Pitch plane contours of Mach number (contour levels of 0.01, 1.0, and 24.4 are shown) for 70° 
sphere-cone with 3 different values of nose bluntness ( RJRb = 0.5, 0.8, and 1.1). Also shown are pitch plane 
contours of Mach number for a spherical-section heatshield of the same diameter as the 70° sphere-cone. The 
results are from computations for a freestream velocity of 5.41 km/s and an angle of attack of -15.7°. 
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It is evident that as the nose gets blunter ( RJRb > 1.1), the heatshield acquires more and more character of a 
spherical section. In other words, the spherical section heatshield can be viewed as the limiting case of a sphere-cone 
heatshield for which there is no frustum! 

Both laminar and turbulent computations are performed for the configurations shown in Fig. 11 using modeling 
assumptions outlined in Table 2 and at the peak heating point along the retrograde/steep trajectory (Table 4). Pitch 
plane Mach number contours of 0.01, 1.0, and 24.4 for the sphere-cone (with different RJRb ratios) and the spherical 
section heatshields are shown in Fig. 12. 

For the sphere-cone geometries, blunting the nose does indeed increase the shock standoff distance, which 
should translate into reduced heating on the apex of the cone, and a reduction in the curvature of the shock at the 
apex of the cone (Figs. 12a through 12c). Furthermore, regardless of the amount of blunting, the sonic line attaches 
at the windward shoulder for all cases, but appears to attach in the vicinity of the sphere-cone junction on the 
leeward side (Figs. 12a through 12c). Therefore, one would expect to see the inflection in pressure on the leeward 
side move outboard towards the shoulder with increasing nose bluntness. For the spherical section, however, the 
sonic line attaches at the windward and leeward shoulders (Fig. 12d). 

Centerline distributions of hot- wall (s = 0.85) heat flux and pressure are shown in Fig. 13. 
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Figure 13. Laminar and turbulent centerline distributions of hot-wall (e =0.85) heat flux and pressure for a 
4.7 m diameter 70° sphere-cone at a flight Mach number of 24.4 and at an angle of attack of -15.7° for three 
different values of nose bluntness - RJRb = 0.5, 0.8, 1.1. Also shown are heat flux and pressure distributions 
for a spherical section at the same freestream conditions and same attitude. Dashed lines represent laminar 
results. 


First, heating of the apex of the cone decreases with increasing bluntness, and moves more toward the heating 
level of the spherical section (Fig. 13a). However, despite the increased bluntness, there is enough of the conical 
frustum left that the laminar heating levels there are still lower than that of the spherical section heatshield. This is 
simply because a conical section has curvature in the circumferential direction but not in the longitudinal one. 
Although we have not attempted to blunt the nose any further than an RJRb value of 1 . 1 , we speculate that for values 
larger than 1.1, heating levels will be closer to that of the spherical section heatshield everywhere. 

Next, as expected, the inflection in pressure on the leeward flank of the cone moves outboard towards the 
shoulder with increasing blunting of the nose (Fig. 13b). For the case of turbulent flow, pressure inflections (there 
are multiple points of inflection) cause rapid increase in turbulent heating. However, since there is little frustum 
length left with increasing nose bluntness, the “recovery” of pressure to laminar levels does not occur soon enough. 
Therefore, turbulent heating towards the shoulder decreases with increasing bluntness. Also, since the pressure 
distributions (and inflections in them) are similar for various nose radii considered, one infers that the behavior is 
driven by the attachment of the sonic line at the sphere-cone junction. 
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All computations performed thus far have considered a fixed angle of attack of -15.7°. Such an angle of attack is 
necessary to provide lift for the 70° sphere-cone (the LID ratio for the baseline 70° sphere-cone at -15.7° angle of 
attack is nearly 0.24). Such a relatively high value of LID is necessary to obtain increased downrange for the entry 
capsule. For larger LID , one could perhaps increase the attitude of the capsule to the point where one meridian 
(longitude) in the pitch plane is perpendicular to the oncoming stream. For the 70° sphere-cone, this angle of attack 
is -20°. At such an angle of attack, the apex of the sphere-cone acts more like a shoulder (of large radius) and one 
would expect the sonic line to attach there. Increasing the angle of attack beyond this point is not particularly 
productive because the shoulder (small radius) and the backshell experience the oncoming flow leading to excessive 
heating in those areas. It is instructive to carry out a sweep in angle of attack from 0° to 6 C - 90° for the same 
freestream velocity and density considered until now. 

E. 70° sphere-cone: Variation in angle of attack 

For the baseline 70° sphere-cone configuration, four angles of attack - a = -5°, -10°, -15.7°, and -20° - are 
considered in this part of the study. For the freestream conditions at the peak heating point of the retrograde/steep 
trajectory (Table 4), laminar and turbulent computations are performed. As before, pitch plane Mach number 
contours (contour levels of 0.01, 1, and 24.4) are presented in Fig. 14. 



(a) a = -5° 



Figure 14. Pitch plane contours of Mach number (contour levels of 0.01, 1.0, and 24.4 are shown) for the 
baseline 70° sphere-cone at a freestream velocity of 5.41 km/s and four angles of attack - a = -5°, -10°, -15.7°, 
and -20°. 


First, for the lowest angle of attack (-5°), the stagnation point is on the spherical nose cap (Fig. 14a), and for the 
next angle of attack (-10°) is still on the sphere, but very close to the junction of the sphere and cone (Fig. 14b). For 
the other two angles of attack -a — -15.7° and -20° - the stagnation points lie on the conical frustum (Figs. 14c and 
14d). 

Next, the bow shock has the most curvature at the apex for the lowest angle of attack (Fig. 14a), with the 
curvature gradually decreasing with increasing angle of attack. Furthermore, the shock standoff distance 
progressively increases with increasing angle of attack. One would then expect heating at the apex to decrease with 
increasing angle of attack. For all four angles of attack, the sonic line attaches on the shoulder on the windward side, 
but close to the sphere-cone junction on the leeward side (Figs. 14a through 14d). 

Centerline distributions of hot- wall (s = 0.85) heat flux and pressure are shown in Fig. 15. 
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Figure 15. Centerline distributions of hot-wall (£=0.85) heat flux and pressure for a 4.7 m diameter 70° 
sphere-cone at a flight Mach number of 24.4 and at angles of attack of -5°, -10°, -15.7°, and -20° for one value 
of nose bluntness - RJRb = 0.5. Dashed lines represent laminar results. 

As expected, the heating levels on the apex of the sphere-cone heatshield decrease with increasing angle of 
attack - the sphere-cone presents an increasingly blunt body to the flow with increasing angle of attack. At the same 
time, much of the leeward side is moving away from the flow, i.e., the shock standoff distance increases on the 
leeside with increasing angle of attack (Figs. 15a through 15d) thereby reducing the heating. At the lowest angle of 
attack (-5°), turbulent heating of the windward flank is substantial, and in increasing the angle of attack ( a ) from -5° 
to -10° there is about a 50% reduction in windward turbulent heating. Beyond this angle of attack, any decreases in 
windward turbulent heating are marginal because the stagnation point moves from the spherical nose cap to the 
conical frustum (a loss in one principal radius of curvature). 

An examination of pressure distributions (Fig. 15b) shows that the differences between laminar and turbulent 
pressure levels increase with increasing angle of attack, although from the windside up to the sphere-cone junction 
on the leeward side and little beyond that, pressures obtained from laminar and turbulent flow solutions are 
indistinguishable from each other. Beyond the second point of inflection, there is divergence between pressures 
predicted using the laminar model and pressures predicted using the Baldwin-Lomax turbulence model. The most 
likely cause of this difference is the modification to the A + term in the Baldwin-Lomax. This modification perhaps 
acts like a feedback loop especially since A + depends on the turbulent eddy viscosity ji t . Furthermore, the amplitudes 
of pressure inflection increase with increasing angle of attack, especially at the sphere-cone junction, even for 
laminar flow. 

A complete investigation of the 70° sphere-cone requires consideration of a range of Mach numbers, dynamic 
pressures, and angles of attack. Since space here is limited, an attempt is made to study the impact of freestream 
dynamic pressure alone on the turbulent heating environment. 

F. 70° sphere-cone: Variation in freestream conditions 

Three dynamic pressures are considered in this part of the study. For the freestream conditions at the peak 
heating points of three trajectories (Table 4), the dynamic pressure varies by as much as 60% between the peak 
heating point along the retrograde/steep trajectory and the peak heating point along the prograde/shallow trajectory. 
Much of this variation is due to differences in freestream density for the three cases; the freestream velocities differ 
by no more than 6%. Laminar and turbulent computations are performed for each of three dynamic pressure 
conditions. Centerline distributions of hot-wall ( £ = 0.85) heat flux and pressure are shown in Fig. 16. 
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Figure 16. Centerline distributions of hot-wall (£=0.85) heat flux and pressure for a 4.7 m diameter 70° 
sphere-cone at an angle of attack of -15.7° and three freestream dynamic pressure values - 14.8 kPa, 11.65 
kPa, and 9.16 kPa. Dashed lines represent laminar results. 

Regardless of the freestream dynamic pressure, the spatial trends in heating (laminar or turbulent) are identical - 
the leeward side experiences a rapid rise in turbulent heating (Fig. 16a). The magnitude of heating depends on the 
freestream dynamic pressure - the higher the freestream dynamic pressure, the higher the heating - which is hardly 
surprising. 

The pressure distributions (Fig. 16b) are also similar across the three freestream dynamic pressures considered. 
The pressures from laminar and turbulent computations are identical up to the sphere-cone junction on the leeward 
side. Beyond this, on the leeward flank, the pressures from turbulent predictions show multiple inflections, with the 
amplitudes of these inflections increasing with increasing freestream dynamic pressure. 

In all of the preceding discussion, the focus has been on the 70° sphere-cone shape. However, results for the 
alternate baseline shape - the spherical section heatshield - have been shown as well. The spherical section shape 
clearly has much superior aerothermal performance compared to the 70° sphere-cone at the same angle of attack. 
Firstly, there are no inflections in surface pressure because there are no curvature discontinuities over the acreage 
except at the junction of the dish section and the shoulder. Secondly, the streamwise pressure gradient is always 
favorable, which means transition to turbulence is less likely. Thirdly, even if transition does occur on the acreage, 
the turbulent heating levels are much lower than those of the 70° sphere-cone. Therefore, it would appear that the 
spherical section heatshield is a natural choice to replace the 70° sphere-cone heatshield for a Mars entry capsule. 
However, the technology readiness level of the spherical section heatshield for flight through Mars atmosphere is 
low, even though it has flight heritage in Earth atmosphere. Furthermore, the aerodynamic trim characteristics for 
the spherical section heatshield are likely to require more ballast mass than the 70° sphere-cone for a lifting entry 
(Viking and MSL), and compared to a 70° sphere-cone, there is significant loss of packing volume for the spherical 
section heatshield. With the development of systems analysis frameworks such as COBRA (Co-Optimization of 
Blunt-body Re-entry Analysis), 23 it is possible that the spherical section heatshield parameters ( R n and R s ) can be 
optimized to meet multiple objectives such as maximizing the drag area while keeping maximum heat flux to a 
minimum and reducing the amount of ballast required. However, the problem of elevating the technology readiness 
level of a spherical section heatshield will remain, and investments will have to be made to characterize its 
aerodynamics and aerothermodynamics. Since this is true of any shape, regardless of its non-Mars flight heritage, 
the present study explores several other options. 

We seek geometries, other than the 70° sphere-cone or the spherical section, that can address one or more of the 
issues encountered during the detailed study of the flight heritage shapes. Specifically, we seek geometries that do 
not have any curvature discontinuities and develop a favorable pressure gradient on the leeward side, although some 
may at best minimize abrupt changes in curvature as seen in traditional shapes (Fig. 3). 
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IV. Geometry Modifications 

Two approaches have been taken in the present work to construct geometries with no curvature discontinuities. 
In the first approach, the power of a general-purpose nonlinear constrained optimization package is exploited, while 
still retaining the idea of a multi-segment/composite analytic generatrix. In the second approach, the use of non- 
uniform rational B-splines, or NURBS, which form the basis of most modern CAD software, is explored. The 
objective behind the use of NURBS is to replace the composite generatrix with a single segment (if possible). Based 
on the lessons learned from the two approaches, an attempt is also made to combine the best features of the two. 

A. Configurations based on nonlinear constrained optimization 

In the initial form of this approach, the heatshield geometry is still composed of three segments (like the 
traditional sphere-cone). The generatrix is composed of two circle segments corresponding to the nose cap and 
shoulder. However, instead of a straight line that is tangent to the two circles, an elliptical one is used. The nonlinear 
optimization package is then used with the constraints that the elliptical segment satisfies both first and second 
derivative continuity at either end. Details of this procedure are given in the Appendix, where it is seen that second 
derivative discontinuities may be minimized rather than forced to zero. 

In the present study three candidate configurations, named Z, Z2, and Z3, have been constructed using nonlinear 
optimization. The generatrices and curvature distributions for the three configurations are shown in Fig. 17 along 
with the generatrices of the 70° sphere-cone and the spherical section. The closed symbols shown in the figure are 
the lower and upper knot points, z.e., the points at which the elliptical segment is blended smoothly into the circular 
segments for the nose and shoulder, respectively. The coordinates of the knot points are given in Table 6. All 
generatrices are lined up at the plane of maximum diameter so that one can gauge any loss/gain of volume compared 
to the 70° sphere-cone (the first baseline shape) and the spherical section (the alternate baseline shape). Visual 
inspection of the generatrices in Fig. 17 shows Configuration Z to be very similar to the 70° sphere-cone, but convex 
everywhere. Configuration Z2 is an attempt to blunt the nose (the lower knot point is moved outboard and aft) and 
improve shoulder blending. Configuration Z3, developed after examining the results of Configurations Z and Z2, is 
developed as a middle ground between Z and Z2 by suitable adjustment of various constraints. 
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(a) Generatrices (b) Curvature distributions 

Figure 17. Generatrices and curvature distributions for three alternate 4.7 m diameter configurations - Z, 
Z2, and Z3 - with improved continuity in curvature at the knot points (closed symbols). The configurations 
are a blend of a circle, ellipse, and circle respresenting the nose cap, the “frustum”, and the shoulder, 
respectively. Also shown are generatrices for the baseline shapes - the 70° sphere-cone and the spherical 
section heatshield. 


20 

American Institute of Aeronautics and Astronautics 


Table 6. Lower and upper knot point coordinates for Configurations Z, Z2, Z3 


Configuration 

Lower Knot Point 

x, m r, m 

Upper Knot Point 

x, m r, m 

Z 

0.017090 

0.200000 

0.752925 

2.275943 

Z2 

0.010156 

0.229284 

0.680203 

2.306438 

Z3 

0.034758 

0.360423 

0.711792 

2.285000 


The three-dimensional surfaces swept by rotating the five generatrices shown in Fig. 17 around the vehicle roll 
axis (the x axis in the present study) are shown in profile view in Fig. 18. Configurations Z and Z3 are similar to the 
70° sphere-cone, while Configuration Z2 has more of a spherical-section character. 



Figure 18. Three-dimensional profiles of a 70° sphere-cone, Configurations Z, Z2, and Z3 (all obtained 
through nonlinear optimization with constraints), and a spherical section heatshield. 

Laminar and turbulent flow computations are performed at the peak heating point of the retrograde/steep entry 
trajectory (Table 4) with all the modeling assumptions given in Table 2. Pitch plane contours of Mach number 
(contour levels of 0.01, 1.0, and 24.4) for each of the three configurations - Z, Z2, and Z3 - are shown along with 
those for a 70° sphere-cone and a spherical section in Fig. 19. 
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Figure 19. Pitch plane contours of Mach number (contour levels of 0.01, 1.0, and 24.4 are shown) for a 70° 
sphere-cone, a spherical section, and three configurations (Z, Z2, and Z3) obtained using nonlinear 
optimization with constraints on continuity in first and second derivative between segments. Results are from 
computations performed at an angle of attack of -15.7° and a freestream velocity of 5.41 km/s. 


It is perhaps instructive to convert the pitch plane images shown in Fig. 19 into three-dimensional 
representations. The sonic surface for each of the five configurations is shown in Fig. 20. The sonic bubbles for the 
70° sphere-cone and Configurations Z and Z3 are similar, while those for Configuration Z2 and the spherical section 
are similar. From the profiles shown in Fig. 19, one can see that Configurations Z and Z3 have a “cone like” 
character, and Configuration Z2 is closer to that of a spherical section. 







(e) Spherical section 


(a) 70° sphere- (b) Config. Z (c) Config. Z2 (d) Config. Z3 

cone 

Figure 20. Three-dimensional views of the sonic surface (M= 1 surface) for 70° sphere-cone, a spherical 
section, and three configurations (Z, Z2, and Z3) based on nonlinear optimization with constraints. Results are 
from computations performed at a freestream velocity of 5.41 km/s and an angle of attack of -15.7°. The solid 
dark line at the apex of the 70° sphere-cone represents the junction between the spherical nose cap and the 
conical frustum. 
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Configuration Z differs little from the 70° sphere-cone, especially at the shoulder. However, the modest amount 
of curvature added to the dish portion of the heatshield is sufficient to move the bow shock outward and prevent the 
sonic line from attaching close to the apex of the heatshield. Therefore, one would expect the heating on the leeward 
flank of Configuration Z to be lower than that of the 70° sphere-cone, and perhaps no adverse streamwise pressure 
gradient in the immediate vicinity of the apex towards the leeward side. Since the apex of Configuration Z3 is a little 
blunter than that of Configuration Z, a further reduction in turbulent heating on the leeward side is expected. 
Configuration Z2 shows characteristics of the spherical section - the dish portion has more curvature than either 
Configuration Z or Z3, and should have the least heating of the three configurations. Furthermore, since 
Configuration Z2 has a more gradual blend of the dish section with the shoulder, heating at the shoulder is expected 
to be much less than the others. These expectations are borne out by the centerline distributions of hot- wall heat flux 
and pressure shown in Fig. 21. 




Figure 21. Centerline distributions of hot-wall (£=0.85) heat flux and pressure for 4.7 m diameter 
configurations (Z, Z2, and Z3) developed using a nonlinear optimization package. These configurations are 
based on smooth blending between three segments - a circle, an ellipse, and a circle - corresponding to the 
nose cap, the dish, and the shoulder of the heatshield. Also shown are results for the 70° sphere-cone and the 
spherical section heatshields. Dashed lines represent laminar results. 


Due to the longitudinal curvature of the dish portion of the heatshield, and consequent increase in shock standoff 
distance, the heating of the apex decreases for Configurations Z2 and Z3 (Fig. 21a). However, since Configuration Z 
is very similar to the 70° sphere-cone, the decrease in heating at the apex is smaller. But, the increased shock 
standoff distance on the leeward side helps reduce the turbulent heating level. 

The pressure distributions (Fig. 21b) for the three new configurations are interesting. Unlike the 70° sphere-cone 
case, the new blended shapes have smooth pressure distributions, and more importantly, these configurations have 
favorable streamwise pressure gradients over the entire wetted acreage. Favorable streamwise pressure gradients 
help in delaying transition, although this cannot be inferred in the present work because of the assumption of a fully 
developed turbulent flow from the stagnation point. 

Although the heating levels are the lowest for Configuration Z, and the streamwise pressure gradient is 
favorable, the pressure levels are lower than those of the spherical section heatshield. This implies that 
Configuration Z3 will have the lowest drag as well. Although the focus of the present paper is only on turbulent 
aeroheating environments and not aerodynamics, we note that the trim angle of attack for a target L/D could be 
different for the different configurations. 
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B. Configurations based on NURBS 

As mentioned earlier, a second approach towards construction of surfaces without large curvature discontinuities 
is based on non-uniform rational B-splines or NURBS, which are the basis of most modem CAD tools. The concept 
of a NURBS curve is briefly explained here, and more details can be found in the work of Peigl. 24 

Consider the conic segment (planar curve) shown in Fig. 22. The tangents at either end point (A and T) of the 
segment intersect at the point V, and M is the mid-point of the line segment connecting the points A and T. Here A 
can be thought of as the apex of the heatshield, and T the point of tangency of the dish and shoulder toms. 


V 



Figure 22. Construction of an elliptical segment using a non-uniform rational B-spline (NURBS) curve. 

If d(X, Y) represents the distance between two points X and Y in the plane, the parameter p of the NURBS conic 
is defined as: 


<7(1, M) 
d(V,M) 


(7) 


The parameter p is exactly 0.5 for a parabola. For0<p<0.5 one obtains an ellipse, and for 0.5 < p < 1 a 
hyperbola. We note here that the parameter p, which does not seem to have a name, is not the same as the 
eccentricity [e = V(l-Z? 2 /a 2 ), where a and b are, respectively, the semi-major and semi-minor axes] of an ellipse, even 
if the two parameters differ by exactly a factor of 2 for a parabola (i.e., e = 1 and p = 0.5 for a parabola). This aspect 
is evident from Eq. 7 - if p 0, that is d( l,M) 0, then e f 0, because if e were zero, one would have a perfect 
circle. In fact, by the definition of p given in Eq. 7, p ~ 0 is nearly a straight line. We note that it is possible to link 
the NURBS curve to a canonical conic via numerical means, i.e., knowing the coordinates (x,r) of 5 points on the 
NURBS curve, the canonical equation for a general conic can be solved numerically for the various coefficients. 

CFD mesh generation software programs, such as Gridgen , 25 also provide the capability to develop simple CAD 
geometries using NURBS, and this capability is exercised in the present work. The intent behind the use of NURBS 
is to move away from circle-line-circle (or even circle-ellipse-circle) composite generatrices, and replace them with 
a single smooth curve, if possible. In the present study, the shoulder torus is retained, but the remainder of the 
heatshield is replaced by a NURBS curve. 

Two new configurations - Configurations X and Y - are constructed in Gridgen using single elliptical segments 
from the apex of the heatshield to a point on the circular shoulder. The generatrices and curvature distributions for 
these two configurations are shown in Fig. 23 and compared with the 70° sphere-cone and spherical section 
heatshield shapes. It should be noted that neither Configuration X nor Configuration Y has a spherical nose , and the 
packing volumes offered by the various shapes are different. In the present work, for Configuration X, a value of 
0.25 was chosen for p, and for Configuration Y, the value of p was 0.3. The angle of the tangent to the shoulder 
torus was different for the two configurations - 70° for Configuration X, and 66° (obtained from Eq. 2b with the 
spherical section parameters - Rt = 2.358, R n = 5.658, R s = 0.134) for Configuration Y. 
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Axial coordinate (x), m 

(a) Generatrices (b) Curvature distributions 

Figure 23. Generatrices and curvature distributions for two alternate 4.7 m diameter configurations - X 
and Y - based on NURBS curves. Both configurations are a blend of a NURBS ellipse for the dish with a 
circle (representing the shoulder). For Configuration X, the value of p is 0.25, while it is 0.3 for Configuration 
Y. The angle at which Configuration X meets the shoulder torus is 70°, and roughly 66° for Configuration Y 
(based on the parameters, R b , R n , and /? 5 , of the spherical section and Eq. 2b). 

Distributions of curvature (Eq. 3) for the four shapes of Fig. 23a are shown in Fig. 23b. Curvature is a constant for 
sections based on circles/spheres, and curvature is zero for the flat frustum section of the 70° sphere-cone (the 
principal radius of curvature is in the circumferential direction). Although the new configurations have continuous 
curvature distributions except where the elliptical segments meet the circular shoulder, there are cusps at the apex of 
each of these geometries where the larger curvature will lead to increased heating at the apex, if the angle of attack 
is zero. Figure 24 shows a magnified view of the apex of a geometry constructed by rotating a tilted and displaced 
ellipse - the red curve is the actual generatrix, and the black curves is its mirror image about the x axis. Although the 
radial or r axis is tangential to the elliptical segment at the apex, the tilt and eccenticity of the elliptical segment 
control the rate of change of the second derivative, hence change in curvature. 

Their smaller radius of curvature at the apex would pose a problem if capsules based on Configurations X or Y 
were to enter the atmosphere along a ballistic trajectory - the stagnation point would then be on the apex. However, 
for lifting entries this is not a particular problem because the stagnation point moves away from the apex on to the 
acreage, which has much less curvature. 
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Figure 24. Magnified view of the apex of an axisymmetric configuration (Configuration X) built using a 
displaced tilted ellipse (the red curve is the upper half and the black is its mirror image in the lower half). 

The three-dimensional surfaces swept by rotating the four generatrices shown in Fig. 23a around the vehicle roll 
axis (the x axis in the present study) are shown in profile view in Fig. 25. At the junction of the dish and the shoulder 
torus, the 70° sphere-cone and Configuration X have the same slope (70°), while the spherical section and 
Configuration Y have the same slope (66°). 



Figure 25. Three-dimensional profiles of a 70° sphere-cone, Configurations X and Y (both NURBS-based 
surfaces), and a spherical section heatshield. Configuration X and the 70° sphere-cone have the same slope at 
the dish-shoulder interface, while Configuration Y and the spherical section have the same slope at the dish- 
shoulder interface. 

Laminar and turbulent flow computations are performed at the peak heating point of the retrograde/steep entry 
trajectory (Table 4) with all the modeling assumptions given in Table 2. Pitch plane contours of Mach number 
(contour levels of 0.01, 1.0, and 24.4) for each of the configurations - X and Y - are shown along with those for a 
70° sphere-cone and a spherical section in Fig. 26. 
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Figure 26. Pitch plane contours of Mach number (contour levels of 0.01, 1.0, and 24.4 are shown) for 
sphere-cone, a spherical section, and two configurations (X and Y) based on NURBS ellipses. Results 
from computations performed at a freestream velocity of 5.41 km/s and an angle of attack of -15.7°. 


70° 

are 


Neither Configuration X nor Configuration Y has a very curved bow shock at the apex, unlike the 70° sphere- 
cone, and despite the relatively sharp apices. The shock standoff distances for either configuration are larger and one 
expects decreased surface heating. For either configuration, the sonic line attaches closer to the leeward shoulder. 
Consequently, the wall-bounded shear layer will have a subsonic Mach number at the edge, and hence, reduced 
heating on the leeward acreage as well. 

The pitch plane images shown in Fig. 26 are converted into three-dimensional representations. The sonic surface 
for each of the four configurations is shown in Fig. 27. It is interesting to note that while Configuration X and the 
70° sphere-cone have the same slope at the dish-shoulder junction, the sonic bubble has a greater lateral extent for 
Configuration X. Also, the shape of the sonic bubble of Configuration Y is similar to that of the 70° sphere-cone, 
even though Configuration Y and the spherical section have the same slope (66°) at the dish-shoulder junction. The 
large lateral extent of the subsonic bubble for Configuration X implies a higher pressure distribution in the region 
compared to the other geometries. Consequently the drag performance of Configuration X is expected to be better. 
In fact, preliminary analysis of Configuration X using the COBRA 23 framework indicated that Configuration X had 
better aerodynamic performance. 



(a) 70° sphere-cone (b) Configuration X (c) Configuration Y (d) Spherical section 

Figure 27. Three-dimensional views of the sonic surface for 70° sphere-cone, a spherical section, and two 
configurations (X and Y) based on NURBS ellipses. Results are from computations performed at a freestream 
velocity of 5.41 km/s and an angle of attack of -15.7°. The solid dark line at the apex of the 70° sphere-cone 
represents the junction between the spherical nose cap and the conical frustum. 


Centerline distributions of surface heat flux and pressure for the two configurations are shown in Fig. 28. 
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Figure 28. Centerline distributions of hot-wall (£=0.85) heat flux and pressure for 4.7 m diameter 
configurations (X and Y) developed using NURBS ellipses. Also shown are results for the 70° sphere-cone and 
the spherical section heatshields. Dashed lines represent laminar results. 


Barring a heating increase at the apex of the two configurations X and Y, laminar heating over the acreage of the 
two configurations is on par with those of both the 70° sphere-cone and the spherical section (Fig. 28a). Further, the 
heating at the apex is about 25% lower than that of the 70° sphere-cone. Turbulent heating over the acreage is on par 
with the spherical section, and about 50% lower than that of the 70° sphere-cone. The improved aerothermal 
performance of these new shapes over that of the 70° sphere-cone can be traced to the centerline pressure 
distributions shown in Fig. 28b. Both Configurations X and Y exhibit an inflection in their pressure distributions at 
the apex. However, the streamwise pressure gradient for both configurations is favorable over the entire leeward 
acreage. Furthermore, the leeward pressure levels are high enough that either configuration offers better drag 
performance than the 70° sphere-cone. For Configuration X, the predicted centerline heat flux and pressure for three 
angles of attack - a = -10.7°, -15.7°, and -20.7° - are shown in Fig. 29. 




Figure 29. Centerline distributions of hot-wall (£=0.85) heat flux and pressure for the 4.7 m diameter 
Configuration X at three angles of attack - a = -10.7°, -15.7°, and -20.7°. Dashed lines represent laminar 
results. 
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Over the relatively large angle of attack range ( Aa = 10°), the aerothermal performance of Configuration X is 
still good. The turbulent heating distribution (Fig. 29a) at the highest angle of attack is reminiscent of the heating 
distributions seen in ground-based experiments (Figs. 5a and 7a). Across all three angles of attack considered, the 
pressure distribution on the leeward acreage shows a favorable pressure gradient, with the inflection in pressure at 
the apex of the heatshield. 


C. Curvature continuity imposed on NURBS-based profiles: Configuration X+ 

Given the advantages offered by Configuration X, an attempt is made to smoothly blend the dish section of the 
heatshield to a shoulder torus by applying the nonlinear constrained optimization capability to a dual-segment 
construction guided by the NURBS result. This should reduce the spike in heating at the leeward shoulder, and at the 
windward shoulder as well. A sequence of optimizations is actually performed by adjusting the constraints to 
approach the shape of Configuration X while imposing curvature continuity where the two ellipses meet (and hence 
everywhere , including at the shoulder). The new blended configuration (called Configuration X+) is compared with 
the baseline (Configuration X) in Fig. 30. 



~v 
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(a) Profiles (b) Curvature distributions 

Figure 30. Profiles of Configuration X (elliptical section with p = 0.25 and tangent to a circular torus), and 
Configuration X+ (a blend of two elliptical segments - one for the dish and the other for the shoulder). 
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(a) Heat flux (b) Pressure 

Figure 31. Comparison of centerline pressure and hot-wall heat flux distributions for Configurations X and 
X+. Results are from flow computations (laminar and turbulent) at a flight speed of 5.4 km/s and -15.7° angle 
of attack. Dashed lines represent laminar results. 
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The centerline distributions of pressure and hot-wall heat flux for Configurations X and X+ at a flight speed of 
5.4 km/s and -15.7° angle of attack are compared in Fig. 31. While the pressures are similar on the windward side 
for both configurations, there is a faster decrease of pressure on the leeward side of Configuration X+, which means 
Configuration X+ will have a lower drag than Configuration X. However, the impact of blending curvature between 
the shoulder and dish section is clearly seen in the heating distributions, especially on the leeward shoulder - 
Configuration X+ has more than 25% reduction in heating on the leeward shoulder, and the distribution is much 
smoother. Although not as substantial, there is reduction on the windward shoulder as well. 

V. Conclusions 

Despite the rather limited focus of this study of the 70° sphere-cone geometry at an angle of attack , the 
following broad conclusions can be drawn: 

1. The discontinuity in curvature at the junction between the spherical nose and the 70° conical frustum directly 
affects surface pressure. Pressure distributions develop an inflection at this junction, and create an adverse 
streamwise pressure gradient locally. Blunting the nose only serves to move the discontinuity outboard 
towards the leeward shoulder. Although results are not shown here in the present work, the inflection in 
pressure can be seen in all sphere-cone geometries, regardless of the cone angle. 

2. The magnitude of the amplitude of pressure inflection at the sphere-cone junction depends directly on the 
proximity of the bow shock to the body, i.e., compression of the shock layer due to thermochemistry draws 
the shock closer to the body. If the bow shock is very close to the body, as it is in flight, the sonic line 
attaches at the sphere-cone junction, thus making the edge Mach number supersonic for the wall-bounded 
shear layer on the leeward flank of the cone. Such behavior is not observed in ground-based tests because the 
flow medium is either a high-Mach number (8-10) ideal gas with a limit of 6 on the density ratio across the 
shock, or a low-Mach number (4-6) real gas, and for either case the shock standoff distance is significantly 
larger than in flight. 

3. For the flight cases, results of simulations assuming laminar flow show that the pressure “recovers” from the 
inflection at the sphere-cone junction to a nearly constant level over the conical flank on the leeward side. 
However, results of simulations performed with a turbulence model (Baldwin-Lomax) show that the pressure 
distributions have inflections beyond the sphere-cone junction, and the pressure distributions for the laminar 
and turbulent cases differ significantly on the leeward flank. The multiple inflections in the pressure 
distribution on the leeward flank of the cone do correlate with changes in the spatial gradients of heating 
towards the leeward shoulder. 

Viable alternate shapes have been developed using two different approaches. In one approach, the concept of a 
multi-segment/composite heatshield is retained, i.e., the heatshield is assumed composed of a spherical nose cap and 
shoulder torus. However, in a departure from the 70° sphere-cone geometry, which joins the nose and shoulder with 
a conical frustrum, these alternate shapes are joined by a “cone-like” surface but with two principal radii of 
curvature (the regular right circular conical frustum has only one principal radius of curvature). The power of 
nonlinear constrained optimization is used to blend the various segments better by reducing (if not eliminating) 
curvature discontinuities between the segments. The configurations developed in this approach seem to fix most of 
the issues with the 70° sphere-cone, and show superior aerothermal performance. 

Viable alternate shapes have also been developed using NURBS (non-uniform rational B-spline) curves, with the 
express intent of moving away from multi-segment/composite heatshields. Configurations developed using this 
approach also fix issues with the 70° sphere-cone and show superior aerothermal performance. Constrained 
optimization was also applied to capture the essentials of a NURBS-based result while replacing the shoulder circle 
with an ellipse to achieve curvature continuity everywhere. 

No attempt has been made, nor a claim, that the newly developed shapes have “optimal” performance when 
considering volume and packaging, although a one-off computation using the Co-Optimization of Blunt-body Re- 
entry Analysis (COBRA) 23 framework developed at NASA Ames Research Center indicated that Configuration X 
was close to a Pareto optimal front. 26 

In closing we note that claims about the new shapes, even if “optimized” using sophisticated 
techniques/methods, are all based on simulations, and that there is no substitute for a focused experimental program 
to understand the aerostability characteristics (especially dynamic stability characteristics) of these shapes. 
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Appendix: Blending Conics Via Constrained Optimization 

While traditional sphere-cone generatrices are constructed algebraically, and families of more refined shapes can 
be developed with CAD and grid generation packages such as Gridgen, a convenient alternative is to exploit the 
power of a general-purpose nonlinear constrained optimization package. If highly converged solutions are obtained, 
any geometric constraints (such as matching of first and second derivatives at blend points) are typically satisfied to 
virtually full machine precision. This is comparable to solving sets of nonlinear equations but easier to implement. 
More importantly, some quantity of interest (an objective function) can also be minimized or maximized in the 
process, so the results are optimal in some sense — invariably with at least one of the constraints active (at a bound). 

The variables being optimized are typically coordinates of blend points and coefficients of curves such as conics. 
Equality constraints (such as for ensuring a blend point lies on some curve, or for matching slopes at blend points) 
are simply coded as expressions involving the variables, or the difference between two such expressions, with the 
associated upper and lower bounds specified as zero. Inequality constraints (such as those on conic coefficients that 
ensure the conic is an ellipse) are similarly coded as expressions, with one of the bounds typically specified as 
slightly above or below zero. 

A good gradient-based optimization method has the option to compute partial derivatives of the objective and 
constraint functions via forward or central differencing, so the application program does not need to provide any 
gradient information. In a reusable software framework, the application specifics are confined to the routines that 
read the variables and constraints by name, read the choice of objective function, evaluate the objective and 
constraint functions, and write results in plottable form. 

Configuration Z Construction Details 

Details of the implementation that led to the Configuration Z generatrix are presented as illustration of the power 
and convenience of the constrained optimization approach. 

Blend points N and S on the nose and shoulder 
circles respectively are permitted to slide along the 
circles that are part of the generatrix of the 
underlying sphere-cone, taken to have a 70° half 
angle in this case. Radial coordinate r is synonymous 
with y. A tilted ellipse is substituted for the conical 
flank, adding convexity and some volume. Initially, 
specifying a flank height (maximum distance from 
the ellipse to the straight line of the sphere-cone) was 
envisaged as a parameter defining a family of shapes. 

However, preserving the slope continuity at N and S 
(let alone curvature) proved to be so restrictive that 
this height was not specified in practice. Instead, 
bounds on the coordinates of N and S were specified 
so as to preserve reasonable portions of the circles, 
and the blending was made as good as possible by 
minimizing the [sum of squares] mismatches in 
second derivatives at N and S. 


The input control file shown below shows how nine optimization variables are specified as coordinates (. x N , y N ), 
(x s , ys ) and the coefficients of the flank ellipse x 2 + 2hxy + by 2 + 2gx + 2fy + c = 0. The conditions for this to be 
an ellipse and not a parabola or hyperbola are b - h 2 > 0 and be - f 2 - h(hc -fg) + g(hf- bg) < 0. The “rho” 
inputs are multipliers on possible contributions to the objective function, and only the mismatches in d 2 y/dx 2 were 
employed for Configuration Z. No linear constraints entered the picture, and the eight nonlinear constraints serve to 
make N lie on both the nose circle and the flank curve (with analogous constraints on S), to make that curve an 
ellipse, and to force first derivative continuity at N and S. Miscellaneous controls specify the underlying sphere-cone 
to a capsule forebody module, followed by namelist-type inputs to the optimization package NPOPT . 27 The two 
variable bounds and seven out of eight constraint bounds that were active at the Configuration Z solution are bolded 
in the solution summary that follows the input control file below. 



x, m 


Figure Al. Configuration Z construction. 
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OPTIMIZE Control File (Circle-Ellipse-Circle Blending) 


Mode # Des Vars # L Cons # NL Cons # Des Cycs # Fn. Evals (max) 
190 8 500 999999 

ETA STEPMX EPSOBJ 

0.1 l.E+20 l.E-13 

Design Variables 


# 

VTYPE 

UP 

LW 

WIDTH 

XCNTR 

XMIN 

XMAX 

V 

VS C ALE 

H 

BL 

BU 

1 

' XN 

0. 

1. 

99. 

99. 

99. 

99. 

0.028 

1.000 

l.E-6 

0.01 

0.08 

2 

' YN 

0 . 

1 . 

99. 

99. 

99. 

99. 

0.250 

1.000 

l.E-6 

0.20 

0.50 

3 

'XS 

0 . 

1 . 

99. 

99. 

99. 

99. 

0.750 

1.000 

l.E-6 

0.65 

0.85 

4 

' YS 

0 . 

1 . 

99. 

99. 

99. 

99. 

2.300 

1.000 

l.E-6 

2.20 

2.35 

5 

'BF 

0 . 

1 . 

99. 

99. 

99. 

99. 

0.185 

1.000 

l.E-6 

0.02 

3.00 

6 

CF 

0 . 

1 . 

99. 

99. 

99. 

99. 

0.005 

1.000 

l.E-6 

.001 

0.20 

7 

' FF 

0 . 

1 . 

99. 

99. 

99. 

99. 

-0.015 

1.000 

l.E-6 

-5.0 

1.00 

8 

'GF 

0 . 

1 . 

99. 

99. 

99. 

99. 

-0.300 

1.000 

l.E-6 

-5.0 

1.00 

9 

'HF 

0 . 

1 . 

99. 

99. 

99. 

99. 

-0.300 

1.000 

l.E-6 

-5.0 

1.00 


Objective function contributions 


TYPE MULTIPLIER I l_OB J I2_0BJ Xl_OBJ X2_OBJ 
RHO_XN 0.00 999 999 999. 999. 
RHO_DYDX 1.00 999 999 999. 999. 
RHO_DY2DX2 1.00 999 999 999. 999. 
Linear Constraints 


J LCTYPE BOUNDLO BOUNDUP TLCON ILCON 
Nonlinear Constraints 


J 

NLCTYPE 

BOUNDLO 

BOUNDUP 

XNLCON 

INLCON 

JNLCON 

SNLCON 

1 

'RESID NN' 

0.000 

0.000 

999 . 

999 

999 

10.0 

2 

'RESID NF' 

0.000 

0.000 

999. 

999 

999 

10.0 

3 

'RESID SF' 

0.000 

0.000 

999. 

999 

999 

10.0 

4 

'RESID SS' 

0.000 

0.000 

999. 

999 

999 

10.0 

5 

' ELLIPSE 1 ' 

0.003 

9999.0 

999. 

999 

999 

10.0 

6 

' ELLIPSE2 ' 

-9999.0 

-0.003 

999. 

999 

999 

10.0 

7 

' DERIV1_N ' 

0.000 

0.000 

999 . 

999 

999 

1.0 

8 

' DERIV1 S' 

0.000 

0.000 

999. 

999 

999 

0.1 

Miscellaneous 

Controls : 

Underlying 

sphere- 

cone definition 



1.178787 ! RADIUS_NOSE 

2.357574 ! RADIUS_BASE 

0.1341402 ! RADIUS_SHOULDER 

70. ! HALF_CONE_ANGLE 

35. ! SKIRT_ANGLE 

0.006675 ! SKIRT_LENGTH 

401 ! # discretization points, approximately uniform 

$NPOPTIONS 

LEVELVER = 1, MAJORPL = 10, MINORPL = 0, MINORIL = 1000, NPRFREQ = 0, 
LEVELDER = 0, DIFFERENCE_INTERVAL = -l.E-6, PENPARAM = 0., STEPLIM = 0.01, 
TOLLIN = l.E-6, TOLNLIN = l.E-9, TOLOPT = l.E-7, 

$END 

Output Summary 


Itns Major Minors Step nCon Feasible Optimal MeritFunction L+U BSwap nS condHz Penalty 

78 50 0 1.0E+00 66 ( 8 . 2E-13 ) ( 3 . OE-15 ) 1 . 7449128E+04 40 1.0E+00 _n r c 

EXIT — optimal solution found 

Problem name NLP 

No. of iterations 78 Objective value 1 . 7449128089E+04 

No. of major iterations 50 Linear objective 0 . 0000000000E+00 

Penalty parameter 0.000E+00 Nonlinear objective 1 . 7449128089E+04 

No. of calls to funobj 880 No. of calls to funcon 880 

Calls with modes 1,2 (known g) 66 Calls with modes 1,2 (known g) 66 

Calls for forward differencing 549 Calls for forward differencing 549 

Calls for central differencing 126 Calls for central differencing 126 

No. of degenerate steps 0 Percentage 0.00 

Norm of x 1.4E+00 Norm of pi 1.6E+05 

Max Primal infeas 15 3.6E-13 Max Dual infeas 7 4.7E-10 

Nonlinear constraint violn 1.9E-12 
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Variable 


State 

Value 

Lower bound 

Upper bound 

Lagr multiplier 

Slack 

variable 

1 

FR 

0 . 1709049E-01 

0 . 1000000E-01 

0 . 8000000E-01 


0 . 7090E-02 

variable 

2 

LL 

0.2000000 

0.2000000 

0.5000000 

61680.84 


variable 

3 

FR 

0.7529250 

0.6500000 

0.8500000 


0 . 9707E-01 

variable 

4 

FR 

2.275943 

2.200000 

2.350000 


0 . 7406E-01 

variable 

5 

FR 

0.2440739 

0 . 2000000E-01 

3.000000 


0.2241 

variable 

6 

LL 

0 . 1000000 E -02 

0 . 1000000 E -02 

0.2000000 

856998.1 


variable 

7 

FR 

-0 . 1195725E-01 

-5.000000 

1.000000 


1.012 

variable 

8 

FR 

-0 . 8743908E-01 

-5.000000 

1.000000 


1.087 

variable 

9 

FR 

-0.4802948 

-5.000000 

1.000000 


1.480 


Nonlin 

constrnt 

State 

Value 

Lower bound 

Upper bound 

Lagr multiplier 

Slack 

nlncon 

1 

LL 

- 0 . 693889 4 E -14 



-3211.126 

0 . 6939E-14 

nlncon 

2 

LL 

- 0 . 2498002 E -14 



-3456.580 

0 . 2498E-14 

nlncon 

3 

LL 

- 0 . 193178 8 E -11 



-8913.383 

0 . 1932E-11 

nlncon 

4 

LL 

- 0 . 1804112 E -14 



-11049.88 

0 . 1804E-14 

nlncon 

5 

FR 

0.1339079 

0 . 3000000E-01 

99990.00 


0.1039 

nlncon 

6 

UL 

- 0 . 3000000 E -01 

-99990.00 

- 0 . 3000000 E -01 

-65221.14 

-0 . 3609E-12 

nlncon 

7 

LL 

- 0 . 188293 8 E -12 



12989.36 

0 . 1883E-12 

nlncon 

8 

LL 

- 0 . 7278622 E -13 



333845.9 

0 . 7279E-13 


DESIGN VARIABLES AT DESIGN 

1 1.7090485686128E-02 
5 2.4407387290302E-01 
9 -4.8029478658902E-01 


ITERATION 65 

2 2.0000000000000E-01 

6 1.0000000000000E-03 


3 7.5292501222252E-01 

7 -1.1957246451232E-02 


4 2.2759428039880E+00 

8 -8.7439076377836E-02 


DYDX_N 

D2YDX2_N 

D3YDX3_N 

DYDX_S 

D2YDX2_S 

D3YDX3_S 

DYDX_FN 

D2YDX2_FN 

D3YDX3_FN 

DYDX_FS 

D2YDX2_FS 

D3YDX3 FS 


5. 8084825715694 E+ 00 
-1. 7369234892112 E+ 02 
1.5133334722849E+04 
2. 3507550099207 E+ 00 
-1.2428438212533E+02 
1. 6692116308434 E+ 04 
5.8084825715695E+00 
-1. 2758225157932 E+ 02 
1.2523566088849E+04 
2. 3507 5500992 15E+ 00 
-4.9832874145213E-01 
7.6808432143066E-01 


Contributions to the objective: 


RHO_XN 
RHO_DYDX 
RHO DY2DX2 


0.000000000000000E+00 
5.652379684678066E-25 
1.744 9 12 808922 693E+04 


NLCON: Nonlinear constraint check: 


# 

NLCTYPE 

INLCON XNLCON 

CURRENT VALUE 

LOWER BOUND 

UPPER BOUND 

VIOLATION? 

1 

RESID_NN 

999 

999.0000 

0.000000 

0.000000 

0.000000 

* * * 

2 

RESID_NF 

999 

999.0000 

0.000000 

0.000000 

0.000000 

* * * 

3 

RESID_SF 

999 

999.0000 

0.000000 

0.000000 

0.000000 

* * * 

4 

RESID_SS 

999 

999.0000 

0.000000 

0.000000 

0.000000 

* * * 

5 

ELLIPSE1 

999 

999.0000 

0.013391 

0.003000 

9999.000000 


6 

ELLIPSE2 

999 

999.0000 

- 0.003000 

-9999.000000 

- 0.003000 

* * * 

7 

DERIV1_N 

999 

999.0000 

0.000000 

0.000000 

0.000000 

* * * 

8 

DERIV1_S 

999 

999.0000 

0.000000 

0.000000 

0.000000 

* * * 


PROBLEM PARAMETERS AT DESIGN ITERATION 65 


Objective function 
Coef. B of flank ellipse 

Coef. C of flank ellipse 

Coef. F of flank ellipse 

Coef. G of flank ellipse 

Coef. H of flank ellipse 


Initial 
2. 0281882612 E+ 06 
1.8500000000E-01 
5.0000000000E-03 
-1.5000000000E-02 
-3.0000000000E-01 
-3.0000000000E-01 


Current 

1.7449128089E+04 

2.4407387290E-01 

1 . 0000000000 E -03 

-1.1957246451E-02 

-8.7439076378E-02 

-4.8029478659E-01 


Total # objective calculations: 881 


Normal termination. 


OPTIMIZE: CPU secs, for this run: 0.02 


Configuration X+ Outline 


The goal of Configuration X+ was to match Configuration X as well as possible while enforcing curvature 
continuity everywhere, mainly to reduce the heating spike at the shoulder. Two tilted ellipses were blended 


33 

American Institute of Aeronautics and Astronautics 



numerically with a variation of the process described for Configuration Z. The only objective function implemented 
was to minimize (or maximize) the x coordinate of the maximum diameter point U on the shoulder. Since the 
circular shoulder of Configuration X is in the forward position of the simple spherical-section, x v was minimized, 
although in practice x v was also an optimization variable that encountered its lower bound of 0.565 m and stayed 
there. 

The nose ellipse had only three coefficients to solve for (since it had to be vertical at the origin), while the outer 
ellipse blending with it at some point T (x T , yr) had the maximum of five coefficients, making 1 1 optimization 
variables total. There were 10 nonlinear constraints: two conditions each for the nose and shoulder conics to be 
ellipses, one for the nose ellipse to pass through T 7 , two for the shoulder ellipse to pass through T and £/, one for it to 
have zero slope at U, and two for matched first and second derivatives at T. 

Since the shoulder was no longer a circle, a heuristic bridge to a notional aft body cone (35° angle) was 
constructed from point U by means of a cubic defined by two points and two slopes. Matching the aft segment’s 
curvature to that of the shoulder ellipse at U was not attempted and not necessarily desirable anyway. 

A sequence of optimizations was actually performed with adjustments to the constraint bounds and/or variable 
starting guesses made in order to nudge the optimized result close to Configuration X (as determined graphically — 
see Fig. A2) while achieving tight convergence to each optimum. 



Figure A2. Comparison of Configuration X+ fit (ellipse-ellipse) with Configuration X 
(ellipse-circle) and the traditional shapes. 

It is hoped that the relative ease of implementation of the constrained optimization is apparent from the 
preceding discussion. The primary advantages of the scheme are: (1) there are no awkward solutions of nonlinear 
equations; (2) the entire input scheme is table driven (seen in the input listing); and (3) it is possible to satisfy 
specified constraints to any desired accuracy. The discussion could prove useful to any future shape optimization 
studies. 
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